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ABSTRACT 



Aims. We investigate the formation of the H i Lyman a and Lyman p lines in an equatorial coronal streamer. Particular attention is paid to 
frequency redistribution for the scattering of the incident radiation. The properties of the spectral lines are studied. 

Methods. The coronal model is given by a global 2.5-D three fluid solar wind model with a particles. The emergent intensities and line profiles 
are calculated from the solution of the statistical equilibrium and radiative transfer equations for an hydrogen atom with 1 1 energy levels under 
non local thermodynamic equilibrium. The formation of the lines results from radiative excitation, collisional excitation, and takes into account 
the coupling with all other transitions between the hydrogen energy levels. 

Results. We present new estimates of the radiative and collisional contributions of the Lyman line intensities within the streamer. It is also 
shown that within the streamer, the full width at half-maximum (FWHM) of the Lyman yS line is a better indicator of the plasma temperature 
than that of Lyman a. These results show that care should be taken when inferring the proton temperature from the Lyman a line profile as 
observed in coronal streamers, e.g. by the Ultraviolet Coronagraph Spectrometer or the Solar Ultraviolet Measurements of Emitted Radiation 
experiments on board the Solar and Heliospheric Observatory. 
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,^ 1. Introduction 



A large amount of our knowledge of the solar extreme ultravi- 
^ olet (EUV) corona comes from the observations of the first two 
Hi Lyman lines at 1215.67 A and 1025.72 A. The Lyman a 
line is indeed the brightest line in this range of radiation. To 
infer the coronal plasma properties from observations it is nec- 
essary to get the line profiles with enough spectral resolution. 
While chromospheric profiles show a self -reversal at line cen- 
tre, the coronal profile of the scattered Lyman a line is close to 
a gaussian profile. 

From the knowledge of the line width one can derive an 
equivalent temperature (ie. including all causes of line broad- 
ening), and with further assumptions, the kinetic (thermal) tem- 
perature of the hydrogen atoms. If the coupling due to charge 
exchange between hydrogen atoms and ions is strong enough, 
then the proton temperature is e qual to the hydrogen tem- 
perature. All en et al ] il998', 2000) concluded from fast solar 
wind mode l s that this coupling is strong up to 3 R©, while 
lOlsen et"al studied slow wind models and found strong 

coupling up to 10 R0 in some conditions. The decoupling 
of hydrogen atoms and protons will occur at different densi- 
ties, depending on which physical assumptions are made in the 
models. In our case, this coupling is strong enough at densi- 
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ties above 10^ cm"^. This cond ition is met wi thin the whole 
streamer (see for further details IlI et aljEoo3) . Thus we can 
obtain a reliable estimate of the proton temperature from the 
full width at half maximum (FW HM) of the Lyman line pro- 
files (see also lMarsch et alJl999t) . 

In this paper we show that in the case of coronal stream- 
ers, radiative transfer calculations in non local thermodynamic 
equilibrium (NLTE) are useful to predict the properties of the 
Lyman a and Lyman (3 lines. The temperature diagnostic from 
these two lines is examined. We describe the streamer model 
in section|3and the radiative transfer calculations in section|3 
Results are discussed in section|3 

2. The coronal model 

The coronal model is obtained from a global three fluid solar 
wind model wi th a particles. It is described in another paper by 
iLi et alJ ( l2006l) . It is a 2.5-D, axially symmetric model, where 
all variables depend on two spatial coordinates, but where 
the three components of vector quantities are retained. In the 
streamer (the closed magnetic field region) no external heat- 
ing is applied. A hot coronal boundary, electron heat flux and 
Coulomb coupling lead to a non-isothermal streamer in which 
all three species (namely electrons, protons, and cc's) have the 
same temperature. The properties of the streamer model that 
we use for this study are shown in Fig. ^ for two heights. One 
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Fig. 1. Streamer models at a distance of 0.05 Rq (left) and 0.88 
Rq (right) above the solar surface. Top panels: proton density 
in cm"^; bottom panels: electron temperature in K. All quanti- 
ties are plotted against the location along the line-of-sight ex- 
pressed in units of the solar radius R0, where the origin is at 
the centre of the LOS. 

has the line-of-sight (LOS) centre situated at a heliocentric dis- 
tance of 1 .05 R0 at the equator (left column) and the other one 
is at a heliocentric distance of L88 R0, slightly below the cusp 
which is located at ~ 2 Rq. The LOS is chosen to be 4 Rq long. 
The top panels in Fig.^present the density variations along the 
LOS of the protons. The electron density is almost the same as 
the proton density due to the small amount of a's. The bottom 
panels of the figure show the temperature variation along the 
LOS. For the present investigation we have made use of 10 dif- 
ferent LOS for a streamer axis along the equator All LOS are 
4 Rq long, but the centre of the LOS is located at an increasing 
distance from the sun. As one goes further away from the sur- 
face, the temperature profiles flatten and the density variations 
become smoother as well as less sharp. 

3. Computation of the emergent radiation 

In order to compute the radiation emitted by the corona 
we use a numerical code which solves the radiative trans- 
fer (RT) equations and the statistical equilibrium (SE) 
equations in NL T E. Th i s code has been descr i bed i n 
nniittehroze etTll dlQQ^: Inouttebroze & T ,abr(^ tOOd: 
Gouttebroze & Heinzell ( 120021) . Here we recall its most rele- 
vant features. The SE is solved for a 10 level -i- continuum hy- 
drogen atom. The RT is solved for the lines and continua in a 
ID pl ane-parallel geometry using a Feautrier method ( Feautrier 
11964 with variable Eddington factors. The boundary condi- 
tions for the RT equations are determined by the radiation com- 
ing from the disk. This incident radiation is allowed to freely 
penetrate the structure. The Lyman a profile from the solar disk 
is taken from OSO-8 observations llGouttebroze et aL 1978). 
For Lym an 13 to Lvman-9 t hey are taken from SUMER obser- 
vations ( Warren et alJl99^ . Thomson scattering and Rayleigh 
scattering are taken into account in the computation of the con- 
tinuous absorption coefficients. 



Fig. 2. LTE departure coefficient for the ground level of hy- 
drogen as a function of heliocentric distance (solid line). The 
'plus' signs give the value of 1 /Ne in arbitrary units (see text). 

It is necessary to perform the calculations in NLTE. 
Because of the importance of the incident radiation emitted by 
the solar disk and the low coronal densities, LTE cannot be 
reached. We define the LTE departure coefficient of the en- 
ergy level / by: 



With this definition, the LTE departure coefficient is be - 1 
for the continuum level c. Here A^, and Nc are the populations 
of the bound level / and continuum c. [^J^^^^ is given by the 
Saha-Boltzmann distribution. Figure|2presents the value of the 
LTE departure coefficients bi of the ground level of hydrogen 
for our 10 LOS. It shows that the ground level is far from LTE, 
and that this tendency increases with height. This can be ex- 
plained by the fall-off' in densities, which makes the incident 
radiation even more predominant in the formation of the hy- 
drogen spectrum at higher altitudes. In fact the variation of bi 
closely follows that of 1 /Ne (inverse of the mean electron den- 
sity) with altitude, where the mean of the density is defined as 

- Jq ne(x)dx/L, with L the total lengfli of the LOS. This 
is illustrated in Fig. |2l where we plotted as 'plus' signs l/Ne 
in arbitrary units normalized so that the values of bi and l/Ng 
at 1.05 Rq are identical. The factor of proportionality between 
bi and l/Ne is dependent on temperature. This result is valid 
as long as the ionization balance of hydrogen is governed by 
coUisional ionization and radiative recombination. The excited 
levels also are far from LTE but to a lesser extent. For instance 
we have b2/bi < 10"^. 

Once the coupled system of SE and RT equations is solved 
we compute the emergent intensity in a particular line from the 
knowledge of the line source function, the source function for 
continuous absorption, and the related absorption coefficients. 
The emergent intensity /v(/i) at an angle 6 between the normal 
to the surface and the LOS such that /i = cos 6 is given by: 

ly(tl)^ f'SMC-'^'^dt/tl, (2) 

Jo 
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where Ty is the optical thickness at frequency v. Of course the 
simpUcity of Eq. hides the fact that the computation of the 
total source function S y{Ty) in NLTE is a non-trivial task. Here 
we use the formulation of the equivalent two-level atom. One 
can then express the hne source function as: 

S'y = s*B* + (1 -s*)Jy. (3) 

In Eq. (|3|l, e* and B* account for all the processes that can 
affect the creation and the destruction of photons in the tran- 
sition at frequency v, while (1 - s*)Jy is the scattering term 
of the source fu nction. Expressi ons for g* and B* are given 
in [Gouttebrozel fl98d) (see also Ijefferiesi 1 19681 Chap. 8.1). 
Through Eq. Q a non-local and non-linear coupling between 
the radiation and the plasma arises. 

The scattering integral in eq. ^ is expressed as: 

1 r°° 

Jy^ — \ R(v',V)Jy.dv'. (4) 

<fiv Jo 

(fy is the normalized absorption profile of the line, and R(v', v) 
is the angle-averaged frequency redistribution function. It gives 
the probability density that a photon absorbed at frequency v' 

is re-emitted at frequency v. 

We follow the treatment of iHumnieJ ll 19621) who studied 
four types of redistribution of radiation scattered from mov- 
ing atoms. The redistribution functions with the greatest signif- 
icance ( Mihalasiil978.) are Hummer's case II (radiation damp- 
ing with coherence in the atom's rest frame), and case III (ra- 
diation and collision damping with complete redistribution in 
the atom's rest frame). A Doppler effect arises when deriving 
expressions for the redistribution functions in the observer's 
frame. The case II redistribution function Rji is well suited to 
describe the scattering o f radia tion in resonance lines jMihalas ) 
Il97 8: Cram & Vardavai ll978h . and we use it for the redistri- 
bution in the Lyman lines. In the Doppler core of the line (2-3 
Doppler widths from line centre) this function is close to com- 
plete redistribution (CR D). Coherency effects are more notice- 
able in the line wings (Mil kev&Mihalaslll973i) . This means 
that in the solar corona the coherency effects of Rjj may be im- 
portant at about L5 A from the line centre of the Lyman lines, 
and thus difficult to detect. The case III redistribution function 
Riii is close to CRD over the entire line profile. It is the pre- 
dominant redistribution mechanism in collisionally-dominated 
regions of the solar atmosphere. For studies of the Lyman a line 
profile in the corona, it is generally assumed that the redistri- 
bution function is well described in the coherent scattering ap- 
proximation for an atom with two shar p energy l evels (zero nat- 
ural line wi dth; see e.g.,'Withbroe et al. 1982b: iLi etaPl 19981: 
lAllen et alJ fr998: Cranmer 1998). This is the case I redistri- 
bution function in IHummed jl962l) . ICranmeii lll998^ studied 
case II and found that it makes a small difference relative to 
case I which would be hardly noticed in the observations of 
the Lyman a profile. This is due to the fact that both case I and 
case II redistribution functions are close to CRD in the Doppler 
core of the Lyman lines, which extend to about 1.5 A from 
line centre at coronal temperatures. The differences between 
case I and case II arise in the far wings of Lyman a where the 
Thomson component of the line (chromospheric Lya radiation 



scattered from electrons in the corona) becomes the major con- 
tributor to the line intensity. It would be interesting to see if this 
effect is similar for Lyman /3. 

In this work we use partial redistribution in frequency 
(PRD) to compute the frequency redistribution function in 
Eq. © for all Lyman lines up to Ly-9, assuming isotropic scat- 
tering in the laboratory frame. Our redistribution function is 
therefore a linear combination of Rji and Rjij. The redistribu- 
tion function Rju is taken to be equal to the complete redistribu- 
tion function given by the product ifyify' as in Mihalas ( 1971). 
Defining the branching ratio y - F, /(F, + F^), where F^ and 
Ff are the radiative and collisional damping constants respec- 
tively, we have (Omont et al. 1972): 

R(v, v') = yRiiA (V, v') + ( 1 - r)^v^v' . (5) 

We compute Fc for each Lyman Une at each position along the 
LOS. It is important to take this spatial variation into account as 
we found that it makes a difference in the width of the emergent 
line profiles if it is neglected. Furthermore, since the Lyman a 
line has extended wings, we use a frequency -dependent colli- 
sional damping coefficient jYelnik et alJl981 ). All other subor- 
dinate lines are treated with the standard CRD approximation 
by imposing y = in eq. (|5}. A comparison between CRD 
and PRD computations shows that CRD alone would be a bad 
approximation for the Lyman lines. However CRD cannot be 
neglected, as we find noticeable differences in the line widths 
between the scattering with Ru only (by forcing 7=1) and 
PRD for Lyman /3. 

4. Results and discussion 

Figure |3|presents the resulting integrated intensities of the first 
two Lyman lines as a function of distance from sun centre. The 
intensities are calculated by summing over the computed line 
profiles. The decrease of the intensity with altitude is more 
pronounced for Lyman /3 than for Lyman a, a first indication 
that they relate to the plasma parameters in different ways. The 
Lyman a intensities comp are well with the corn puted and ob- 
served values presented inlVasau ez et alJ ll20Q3h . We find that 
the variation of the Lya intensity with height is not related to 
the variation of the electron density, while the decrease of Ly/3 
and Ha (not shown on Fig.|3|l intensities follows the decrease 
of nl closely up to ~ 1.5 Rq, consistent with the fact that these 
two lines are mostly formed by collisional excitation in the in- 
ner corona. Then the fall-off in intensity is less rapid than that 
of the square of the electron density, owing to the growing im- 
portance of radiative excitation (see also Fig.|5}- In fact there is 
a coupling between Lyfi and Ha, which means that an Ho' pho- 
ton can be absorbed and subsequently lead to an emission of a 
Ly/3 photon. We obtain a nearly constant ratio between Ly/3 and 
Ha intensities, with /(LyjS)//(Ha) ^ 8. This coupling lea ds to 
lowered coherency effects, as was illustrated by Heinzel e t alJ 
lll987t) . The coherence coefficient y (eq. |5|i is close to 1 for 
Lya and 0.57 for LyyS at the centre of the LOS closer to the 
Sun (height of 1.05 Rn). We have comp uted the parameter 
A - (Aji/Pj) X y as in H einzel et 

ID CHI 

- although we use 

a slightly different definition for y, where Aji and Pj are the 
spontaneous emission coefficient in the j — » i transition and 
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Fig. 3. Integrated intensities in photons s ' cm"^ sr ' of 
Lyman a (triangles) and Lyman /3 (squares) as a function of 
heliocentric distance. 

the total depopulation rate of level j, respectively. Close to the 
Sun, in the streamer base, our value of A is around 0.99 for Lya 
and 0.3 for Ly/3 (at the centre of LOS). This confirms that the 
coherency effects in Ly/3 are less important than in Lya. 

Now we turn to the study of the Lyman line profiles. We 
fit each computed profile with a gaussian profile obtained from 
a non-linear least squares fit. We find that the Lya wings are 
broader than our gaussian fits at all altitudes. The line centre is 
well reproduced by the gaussian fits at altitudes above ~1.6 Rq. 
The Ly/3 line deviates from the gaussian fits as the altitude in- 
creases. The line centre is always well reproduced by the gaus- 
sian fit, but the wings get broader with height. From our gaus- 
sian fits we obtain the full width at half maximum (FWHM) of 
the computed line profiles. As we go further up in the corona, 
the line widths decrease. This reflects the decrease in temper- 
ature. However, care should be taken when inferring the tem- 
perature from the line width. As the profile is not exactly gaus- 
sian, this introduces an error in the temperature derivation (for 
a more thorough discussion see, e.g., Allen et al. 1998). In this 
respect it seems more reliable to exploit the diagnostic possi- 
bilities of the Lyman /? line which is closer to a gaussian than 
the Lyman a line. 

To infer the temperature of the neutrals from the line 
widths, one can relate the FWHM and the temperature Th with: 

^2 



m 
2k 



FWHM' 



4/l2ln2 



(6) 



We arbitrarily chose to use ^ - 20 km s ' for the non- thermal 
motions in the calculations of the line profiles at all altitudes. 
While it might not be accurate, the exact value for ^ is not im- 
portant for our discussion on the temperatures infeiTed from 
the Lyman lines, as we assume that they both have the same 
non-thermal broadening. The resulting temperatures derived 
from the width of the Lyman a line and the Lyman /3 line are 
plotted in Fig. I^as a function of height of the LOS, together 
with the mean temperature derived from the model input and 
the temperature at the centre of the LOS. The mean temper- 

rM 

ature is defined as Tmean — Jy T(m)dm/M, with m the col- 



Fig. 4. Temperature derived from the width of Lyman a (trian- 
gles), Lyman /? (squares), temperature at centre of LOS (thick 
solid line), and mean temperature from the model input (thick 
dotted line), as a function of heliocentric distance. 

umn mass along the LOS, and M the total column mass of the 
LOS. Inspection of figure|4]clearly shows that the width of the 
Lyman /3 line is indeed a good indicator of the plasma tempera- 
ture from 1.05 R0 up to ~2 Rq. Above this height, the temper- 
ature derived from the Lyman /3 FWHM is slightly lower than 
the mean plasma temperature T^iean and than the central tem- 
perature, an indication that the plasma conditions have signifi- 
cantly changed. Indeed, from the streamer model the magnetic 
cusp is located at about 2 Rq. Within that distance, outflow 
velocities are very small (a few km s"'). Above the cusp, the 
protons reach a velocity of about 100 km s"' at 3 Rq. 

Close to the Sun, the temperature derived from the LyyS 
FWHM is in very good agreement with the temperature at 
the centre of the LOS. However, higher in the streamer, the 
temperature derived from the Ly/3 line width is in better 
agreement with the mean plasma temperature than with the 
temperature at the centre of the LOS. Due to the temperature 
and hydrogen density variations along the LOS (see Fig. 
the contribution of coUisional excitation in the formation of 
the LyyS line is more concentrated at the centre of the LOS 
when close to the Sun, and more smoothly distributed along 
the line of sight higher in the streamer. We note that when 
neglecting CRD in frequency redistribution in eq. (|5}, the 
temperature derived from the Lyy6 width exactly matches the 
mean temperature of the models up to a height of 1.6 Rq. 
It can also be seen from Fig. 0]that the temperature derived 
from the Lyman a line significantly underestimates the mean 
plasma temperature and the temperature at the centre of 
the LOS. This diff'erence can be as much as 3.8 x 10^ K at 
r - 1.39 Rq for the mean temperature, and 5.0 x 10^^ K at 
r - 1.58 R0 for the central temperature. The difiference in 
behaviour between the two Lyman lines points to the different 
relative contributions of radiative and collisional excitation in 
the formation of the two lines. From eq. (|3} we identify the first 
term of the right-hand side with the collisional component and 
the second term with the radiative component. Figure|5]shows 
their variation with altitude for the Lyman a and Lyman /3 
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Fig. 5. Relative contribution of the radiative (solid line) and col- 
lisional (dotted line) components for Lyman a (triangles) and 
Lyman yS (squares) as a function of heliocentric distance. 



lines. The radiative component is represented with a solid 
line, and the collisional component with a dotted line, while 
triangles stand for Lyman a and squares for Lyman /3. This 
figure shows that, as expected, Lyman /? is mostly formed by 
collisional excitation in the streamer. However it is interesting 
to note that the collisional component of Lyman a is not 
negligible in the streamer base, contributing to nearly 10% 
of the line intensity close to the Sun. We should stress here 
that both the scattered and collisional components will include 
contributions resulting from exchanges with other transitions 
(lines and continua) involving the rest of the atomic states. The 
va lues reported in Fig. Bl can be compared with those given 
bv lRavmond et al.l(l 19971) in the centre of a streamer observed 
by UVCS. These authors obtained a collisional contribution 
of 1.1% for Lyman a and 57% for Lyman /3 at log T = 6.2, 
while our values at that temperature (at a heliocentric height of 
1.7 Rq) lead to 3% and 77%, respectively. 

In this paper NLTE radiative transfer calculations are per- 
formed to compute the properties of the Lyman lines of hy- 
drogen in the solar corona. It is shown that the width of the 
Lyman /3 line is a better indicator of the plasma temperature 
than the width of the Lyman a line, especially within the 
streamer. It is due to the formation mechanisms of the lines, 
and the coupling of Ly^S with Ho-. This work has been done 
using the approximation of isotropic frequency redistribution 
in the laboratory frame. We do not expect that the inclusion 
of angle-dependent redistribution functions in our calculations 
would change the main conclusion of this Letter, namely that 
the width of Lyman /? is a better proxy for the plasma tem- 
perature in the streamer. It is known that the consideration 
of angular redistribution with dipole scattering has the effect 
of narrowing the line profile compared to isotropic redis tribu- 
tion, due to non-90° scatt e ring of photons (e.g. Withbroe et alJ 
Il982at lAllen et al.1 [19981 ICranm er 1998). Therefore we be- 
lieve that the discrepancy found here between temperatures de- 
rived from Lyo- line profiles and model temperatures could be 



even greater by using angle-dependent redistribution functions. 
Furthermore. ICram & Vardavasllll978ll showed that the effects 
of angle-dependent PRD will be more important when the in- 
cident lines show substantial center-to-limb variations. This is 
not the case of hydrogen Lyman lines, but this is the case for, 
e.g., the O VI lines at 1032 and 1038 A. We will include this ef- 
fect in a future study when we include O vi in our calculations. 
We also obtain new estimates of the radiative and collisional 
contributions of the Lyman line intensities in a non-isothermal 
streamer. These new values may have some importance in the 
derivation of element abundances. Element abundances rela- 
tive to hydrogen can be inferred independently from the ra- 
tio of the resonantly scattered (or collisional) component of a 
spectral line to the res onantly scattered (or collisional) com- 
ponent of, say, H i Ly^S dWithbroe et aljfl982at iRavmond et alJ 
1997). Therefore the relative contributions of the two com- 
ponents of the Lyman lines to their total observed intensities 
must be known with good accuracy. Finally, our results can be 
compared with observations by the SUMER and UVCS spec- 
trometers on SOHO. The radiative transfer calculations can be 
enhanced by including other effects such as Doppler dimming 
to improve the modelling in regions of the corona where out- 
flow velocities canno t be ignored. This has been presented in 
iLabrosse et 1111120061) . 

Acknowledgements. The authors are grateful to the referee. Dr. 
P. Heinzel, for his thorough comments that improved the clarity of 
this work, and to P. Gouttebroze, S. Habbal and J.-C. Vial for their 
critical reading of an early version of the manuscript. Support from 
PPARC grant PPA/G/O/2003/00017 is acknowledged. 

References 

Allen, L. A., Habbal, S. R., & Hu, Y. Q. 1998, J. Geophys. Res., 
103,6551 

Allen, L. A., Habbal, S. R., & Li, X. 2000, J. Geophys. Res., 
105,23123 

Cram, L. E. & Vardavas, 1. M. 1978, Sol. Phys., 57, 27 
Cranmer, S. R. 1998, ApJ, 508, 925 
Feautrier, R 1964, C. R. Acad. Sci. Paris, 258, 3189 
Gouttebroze, P. 1980, PhD thesis, Faculte des Sciences de Paris 

7, Universite Paris 7 
Gouttebroze, P & Heinzel, P 2002, A&A, 385, 273 
Gouttebroze, P, Heinzel, P, & Vial, J.-C. 1993, A&AS, 99, 513 
Gouttebroze, P & Labrosse, N. 2000, Sol. Phys., 196, 349 
Gouttebroze, P, Lemaire, P, Vial, J. C, & Artzner, G. 1978, 

ApJ, 225, 655 

Heinzel, P, Gouttebroze, P, & Vial, J.-C. 1987, A&A, 183, 351 

Hummer, D. G. 1962, MNRAS, 125, 21 

Jefferies, J. T. 1968, Spectral Une formation (A Blaisdell 

Book in the Pure and AppUed Sciences, Waltham, Mass.: 

Blaisdell, 1968) 
Labrosse, N., Li, X., & Li, B. 2006, in ESA SP-617: SOHO-17: 

10 Years of SOHO and Beyond 
Li, B., Li, X., & Labrosse, N. 2006, J. Geophys. Res., Ill, in 

press, doi:10.1029/2005JA011303 
Li, X., Habbal, S. R., Kohl, J., & Noci, G. 1998, ApJ, 501, L133 
Marsch, E., Tu, C.-Y, Heinzel, P, Wilhelm, K., & Curdt, W. 

1999, A&A, 347, 676 



6 



N. Labrosse et al.: On the Lyman a and Lyman /3 lines in solar coronal streamers 



Mihalas, D. 1978, Stellar atmospheres /2nd edition/ (San 

Francisco, W. H. Freeman and Co., 1978. 650 p.) 
Milkey, R. W. & Mihalas, D. 1973, ApJ, 185, 709 
Olsen, E. L., Leer, E., & Holzer, T. E. 1994, ApJ, 420, 913 
Omont, A., Smith, E. W., & Cooper, J. 1972, ApJ, 175, 185 
Raymond, J. C, Kohl, J. L., Noci, G., et al. 1997, Sol. Phys., 
175, 645 

Vasquez, A. M., van Ballegooijen, A. A., & Raymond, J. C. 

2003, ApJ, 598, 1361 
Warren, H. R, Mariska, J. T., & Wilhehn, K. 1998, ApJS, 119, 

105 

Withbroe, G. L., Kohl, J. L., Weiser, H., & Munro, R. H. 1982a, 

Space Science Reviews, 33, 17 
Withbroe, G. L., Kohl, J. L., Weiser, H., Noci, G., & Munro, 

R. H. 1982b, ApJ, 254, 361 
Yelnik, J.-B., Burnett, K., Cooper, J., Ballagh, R. J., & 

Voslamber, D. 1981, ApJ, 248, 705 



